Non-MarkofRan effects of a simple nonlinear bath 
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We analyze a model of a nonlinear bath consisting of a single two-level system coupled to a linear 
bath (a classical noise force in the limit considered here). This allows us to study the effects of a 
nonlinear, non-Markofhan bath in a particularly simple situation. We analyze the effects of this 
bath onto the dynamics of a spin by calculating the decay of the equilibrium correlator of the spin's 
2-component. The exact results are compared with those obtained using three commonly used ap- 
proximations: a Markoffian master equation for the spin dynamics, a weak-coupling approximation, 
and the substitution of a linear bath for the original nonlinear bath. 

PACS numbers: 03.65.Yz, 05.40.-a 
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I. INTRODUCTION 

The linear bath of oscillators plays a prominent role in 
discussions of dissipation and decoherenceBa. In the clas- 
sical limit, the force fluctuations derived from that bath 
correspond to a Gaussian random process. Although this 
is a generic case (due to the central limit theorem), there 
are physical situations when non-Gaussian random pro- 
cesses are important. In this report, we examine the 
simplest possible quantum-mechanical bath whose fluc- 
tuations correspond to a classical telegraph noise: a sin- 
gle two-level system subject to a white noise force. The 
effects of this nonlinear bath are analyzed by coupling it, 
in turn, to a spin, whose relaxational dynamics under the 
action of the bath is calculated. 

In the literature, another type of physically tfilevant 
nonlinear bath is usually discussed: the spin bathEra, con- 
sisting of some large number of spins which are coupled 
to the system under consideration. Our model system 
is simpler in that it contains only a single "nonlinear 
element" , the two-level system. Irreversibility is gen- 
erated not by having a larger number of spins but by 
the coupling to the linear bath. Although designed as a 
drastically simplified model system, it may be physicallss 
relevant for situations like charged tunneling systemsB 
in the vicinity of a mesosconic quantum-coherent de- 
vice (e.g. a Cooper-pair boxQ), which lead to electro- 
static potential fluctuations and which are, themselves, 
also subject to dissipation and decoherence by their en- 
vironment. Viewed as a whole, our model consists of 
two coupled two-level systems, one of which is coupled 
to a linear bath. Of course, such systems have been 
studied before, both in the context of the quantum mea- 
surement problem and decoherence of coupled qubit sys- 
tems. In Refs. 1^, ^, ^ and |l^, the model of two spins 
(qubits) coupled to an environment has been analyzed 
in detail. However, we emphasize that our perspective 
and the questions addressed in this article are different 
from these approaches, since we are interested primarily 
in the differences arising from substituting the nonlin- 
ear bath (in the form of the dissipative two-level system) 
by a linear bath (see Fig. [|). This question is relevant. 



since, in many physical situations where the precise na- 
ture of the bath decohering a given system is unknown, it 
is simply treated as a linear bath, with some given corre- 
lation function. It is, therefore, desirable to understand 
in more detail the kind and magnitude of possible errors 
introduced by such an approximation, in cases where the 
coupling cannot be assumed to be very weak. 

The basic strategy is to calculate the equilibrium corre- 
lator of the two-level system exactly (which can be done 
in the limit of infinite temperature) and to compare the 
results to three common approximations. One of those 
involves replacing the nonlinear bath by a linear bath, 
whose correlation function is prescribed to be the same 
as that of the nonlinear bath. The others are a Markof- 
flan master equation and a weak-coupling approximation, 
applied to the dynamics of the spin under the influence 
of the bath. 

The remainder of this work is organized as follows: In 
Sec. H we give the model Hamiltonian. In Sec. pi|, 
the four different approaches are defined by specifying 
the evolution equation for the density matrix in each 
case. Afterwards, we explain (Sec. IV) how the equi- 
librium correlator of the z-component of the spin may 
be obtained by solving these equations. Finally, Sec. 0. 
presents plots showing the numerical results for the spin 
correlator, along with a comparison between the different 
approaches. 



II. THE MODEL 

We consider a two-level system S coupled to another 
two-level system i?, which represents an example of a 
nonlinear dissipative bath, since it is subject to a fluctu- 
ating force F: 
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Jo-^o-f -t- Ao-^ + Fo-f + i/F- (1) 



Here the parameters 65 and A5 serve to define any de- 
sired two-level system S. This system is coupled to B via 
(jf , with the coupling strength between S and B being 
given by J. The oscillations of (jf{t) at the frequency 




will discuss below. The limit of infinite temperature is 
dictated mostly by the desire to have a comparatively 
strong decay of the correlator of <jf{t) (with a decay 
rate on the order of B's transition frequency 2 A) while 
still retaining the validity of a simple Markoffian master 
equation description (for the full system S + B). 

In the following, we will call the exact solution "ap- 
proach 1" , while "approach 2" refers to a simple master 
equation applied to S alone, "approach 3" replaces the 
nonlinear by a linear bath, and "approach 4" is the weak- 
coupling approximation. 
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FIG. 1: (a) Schematic representation of two stochastic pro- 
cesses corresponding to a classical two-level fluctuator, or 
"telegraph noise" (top), and a Gaussian process (bottom), 
yielding the same power spectrum (right), (b) In our model, 
the two-level fluctuator B is coupled to a noise force F and 
therefore represents a (nonlinear and non-Markoffian) bath 
that acts on a system S. (c) The exact master equation de- 
scription ("approach 1" in main text) treats S and B as a 
composite system, subject to F. 



2A are noisy, due to the action of the fluctuating force 
_F, which may stem from a linear bath of harmonic os- 
cillators, whose Hamiltonian is given by Hp. Below we 
will specialize to the limiting case of infinite temperature, 
where F becomes a purely classical noise force. 

The dissipative dynamics of S can be characterized in 
terms of several different quantities. Here we will analyze 
the decay of the equilibrium correlator ^irf (t)(Tf (0)). 

Solving the full model of a system of two interact- 
ing spins coupled to a linear bath at arbitrary temper- 
atures and coupling strengths represents a formidable 
problem in itself. It has been attacked in thp past 
using the Fgynman- Vernon influence functionalO, both 
analytically j-| (in certain limiting cases) as well as 
numericallyt2l. For our purposes, we will be content with 
analyzing a technically simpler special case. We choose 
F to be a classical white noise fluctuating force, 



III. MASTER EQUATION DESCRIPTION OF 
THE RELAXATION 



Master equation for approach 1: The following exact 
master equation description is used for the action of F 
onto the combined system B + S: 



Pssit) = -i[HsB,PsBit)] - jpssit) + -faf psB{t)a^ . 

(3) 

HsB is the Hamiltonian for the system S and B alone. 
Note that, in contrastj-tp the usual master equation, no 
secular approximationtS has been used in deriving this 
equation, which means the resulting decay rate does not 
have to be small when compared to the transition fre- 
quencies of the system S+B. This is possible because the 
bath correlation function is a delta function, which also 
makes the equation exact. We remark further that Eq. 
(^) is solved directly in the basis in which and are 
diagonal. A transformation into the interaction picture 
(as is commonly performed for the usual master equa- 
tion description) would lead to explicitly time-dependent 
terms in this equation. 

Master equation for approach 2: As has been explained 
above, we will use the master equation description not 
only for the action of F onto the combined system S + B 
(see Eq. (^)), but also for the action of F + B onto S 
alone. This constitutes the approximate "approach 2" , 
involving the usual kind of master equation, which is 
valid only for sufficiently weak coupling J, since it is de- 
rived by applying both the Markoff and secular approx- 
imation (see Ref. In the unperturbed eigenbasis of 
system S, it reads: 



{F{t)Fm=j6{t), (2) 

which corresponds to the limit of infinite temperature 
T of the bath (taken such that the overall noise strength 
remains constant). 

Under these circumstances, the dissipative dynamics 
oi S + B under the action of F can be described exactly 
by using a Markoffian master equation. Note that this is, 
of course, unrelated to the validity of a master equation 
description for the action of B onto S alone, which we 



Pskj = ~{Tk+ Tj + Tkj + i{Ak - Aj) (4) 
+i{Ek-Ej))pskj+SkjJ2Psii\^kif'^^ {BB)j^^_E^ . 

Equation (^) describes the relaxation of the reduced 
density matrix ps of system S alone, under the action 
of the coupling Jafaf to the bath. We have introduced 
the abbreviation A = af. 

The Fourier transform of the correlator of _B = Jaf 
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defines the "bath spectrum" : 

(i?(i)s(o)) ^ (mm) = j^(af (t)af (0)) 

1 /■+°° 

(BB)^^— die-*(i?(t)B(0)) . (5) 

J-oo 

It is real and symmetric in the hmit of infinite temper- 
ature considered here, and therefore it is equivalent to a 
classical colored noise force. As will be explained below, 
{BB)^ is found by applying the master equation (||) to 
B alone. The decay rates are defined by 



n 

ffc, = -27rAfcfc% {BB)^ , 
and the energy shifts are given via 



Afc =^|Afcnpy" duj — 



En 



(6) 



(7) 



IV. RELAXATION OF DENSITY MATRIX AND 
CORRELATOR 



Decay of the equilibrium correlator - We want to obtain 
the equilibrium correlator of (xf (i), 

(af Waf (0)) ^ tr [p^^ af ^af (0)] . (11) 

It is convenient to rewrite Eq. (^l]) in terms of the 
projector onto the spin-up state of 5, P = It)^ (TIs = 



(af(t)af(0))=4(P(t)P(0) 



(12) 



Here, we have used p^g^ — \. The correlator of P{t) 
can be found by calculating the probability to find the 
system S in the state "up" at the time if it had been 
"up" at time 0. This has to be averaged over all realiza- 
tions of the random process F{-): 



Here the indices and energies refer to the unperturbed 
eigenstates of the original Hamiltonian of S alone: Hs = 



Approach 3 consists in replacing the nonlinear bath by 
a linear one. If the two-iepil fluctuator B were replaced 
by a harmonic oscillatorEjEil, this procedure of substitut- 
ing a linear bath with an appropriate correlation function 
for the combination of F and B would be exact. Here, 
it is an approximation whose reliability we want to ana- 
lyze by comparison to the exact solution. In our case, the 
fact that the power spectrum (BB)^ is real and symmet- 
ric means that B can be treated as a classical Gaussian 
random process. Therefore, we have to solve 



Psi^) = -■i[Hstoch{t),ps{t)] 



(8) 



with the stochastic time-dependent Hamiltonian: 

Hstochit) = esaf 4- Agaf + B{t)a^ . (9) 

The numerical procedure used for averaging over the re- 
alizations of the process B{-) is described below, in Sec. 

Approach J^: Instead of the Markoff approximation one 
can use a weak-coupling approximation. This keeps the 
full information contained in the correlator (BB)^, at 
the price of introducing a kernel for the master equation 
which is no longer local in time. We find the following 
weak-coupling equation, where only terms up to second 
order in J have been kept: 



Psit) = -i[Hs,ps{t)] 
dT {B{t)B{0)) 



al,ps{t-T) 



(10) 



This equation is conveniently solved by using the Laplace 
transform. 



P(t)P(0) 



-,trB{{r\sUF{t)P^P^s'^uUt)n)s 



F 

(13) 

Here Upit) is the time-evolution operator for S + B 
under the action of a given realization of F(-). This equa- 
tion is valid only because, in our model, the probability 
of finding "spin up" at a certain instant of time is inde- 
pendent of the history of F{-). The expression ( p^ ) is 
nothing but the population psuit) of the state It)^ for a 
time-evolution starting from the initial condition of "spin 

up", P5b(0) = P< 



P{t)P{0)) = -psnm. 



(14) 



Note that psn decays towards 1/2, such that 
(o-f (t)d'f (0)) vanishes for i ^ oo (as it should be). We 
have used the fact that the correlator is symmetric in 
time, since the potentially antisymmetric imaginary part 
vanishes. ps{t) can be calculated by applying the mas- 
ter equation that describes the action of F onto S + B. 
Put differently, Eq. (^ constitutes an example of the 
quantum regression theorem. Using ps, we calculate the 
Fourier transform of the equilibrium correlator of af (t) : 



2tt 
1 

TT 



4-00 



oo 

+ 00 



dte-*(af(t)af(0)) 
dte^^\psiii\t\)~l/2). (15) 



K^^{uj) is real- valued, symmetric and the integral over 
all frequencies gives 1. 

The master equation for the density matrix p (= psb) 
in the four-dimensional Hilbert space of -I- _B represents 
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a system of linear differential equations with constant 
coefficients. The latter are given by a complex-valued 
16 X 16 matrix C that corresponds to the "superopera- 
tor" on the right-hand side of the master equation. The 
solution is the complex vector p, which consists of the 16 
components of the density matrix p : 



-Cp. 



(16) 



The entries of C can be read off directly from Eq. (|3|) 
The formal solution of Eq. ([l6|), 



Pi.t) 



-ct 



(17) 



can be expressed in terms of the right-eigenvectors |/0^^^), 
the left-eigenvectors (p^^-* | and the eigenvalues A'--'-' of C: 
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FIG. 2: The (BB)^ correlation function for different values 
of 7. The parameters are: J = 1, A = 1 



(18) 



C is not necessarily hermitian, so that the A*^-'-' usually 
are complex-valued (with non-negative real parts) and 
the jp^-'-') do not form an orthonormal basis (however, 
^p(*) I — Sij by construction). In order to obtain 

P5ii(^)j we have to perform the trace over B, psii(t) = 
PSBiiiiit) + psBi2i2it). (In psBs'b'sb the indices s,s' 
refer to S, while b, b' refer to B.) We will use the same 
notation for the components of p^-'-', which is a complex 
vector. Then we obtain 



/"OO 

/ die-*psii(t)-E(^nii- 
Jo 



.,o.„,<«. ,,,, 



Takin g th e real part of this expression gives K^^{u!), 
see Eq. (|l5|). 

An analogous formula holds for the master equation 
describing the decay of ps (= p in that case) under the 
action of + i? (approach 2), see Eq. (^). Then, C 
corresponds to the 4x4 matrix whose entries are read 
off from Eq. (^. Therefore, Eq. ( p^ only contains p^i 
in that case, instead of the sum inside the brackets. 

Similarly, we have to obtain the equilibrium correlator 
of af{t), which is needed as input both for the mas- 
ter equation describing the relaxation of S alone (ap- 
proach 2), the numerical sampling of random processes 
(approach 3) and the weak-coupling approximation (ap- 
proach 4). 

This is done by calculating the relaxation of ps (t) un- 
der the action of -F, starting from the initial condition 
/5b(0) = It)^ (TIb and applying the same formulas as 
above (with B instead of S), for the master equation 
adapted to the two-dimensional Hilbert space of B (with 
a 4 X 4 matrix C). (af {t)(Tf (0)) undergoes damped os- 
cillations. Its Fourier transform 



J2 



TT (Cj2 -4A2)2 +4^272 



(20) 



consists of broad peaks of width 7 (for 72 < 4A2), 
which is proportional to the strength of the noise force 
F and may be comparable to the transition frequency 
2 A itself (see Fig. ||). Thus, B indeed represents a noisy 
two-level fluctuator, which acts onto 5 as a nonlinear 
(non-Gaussian) and non-Markoffian (colored) bath. 

In order to solve the equations (10) of approach 4 
(weak-coupling approximation), we need the Laplace 
transform of the equilibrium correlator of the bath _B, 
which is connected to the Fourier transform in the usual 
way: 



Cbb{s) = / dte-'' {B{t)B{Q)) 
Jo 



(21) 



duj- 



(BB) 
s + iuj 



s + 27 



s2 + 27s + 4A2 



Using the Laplace transform, the system of differential 
equations becomes a system of linear algebraic equations, 
which can be solved by matrix inversion. All the re- 
sults can be obtained analytically. However, here we only 
present the comparatively brief expression for the special 
case of es = 0: 



i^f,(c.) = -3?e{ — 



s + ACbb{s) 



4Cbb(s)+4A2 ' 



>}■ (22) 



NUMERICAL RESULTS 



The following steps have been performed for calculat- 
ing the correlator K^^iuj) of (jf (t): 

Approach 1 ("exact description"): The entries of the 
matrix C are obtained from Eq. (H) . The eigenvalues and 
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eigenvectors of C are calculated numerically and used to 
get K^^ioj) according to Eqs. (19) and (|l5|). 

Approach 2 ("simple master equation for S"): First the 
action of F onto B is considered, to obtain the correlation 
function (BB)^. This result is given in Eq. (pO|). It is 
used to set up the master equation describing the action 
of F+B onto S, Eq. (|). Its coefficients define a 4x4 "C- 
matrix" , which is diagonalized. The results are inserted 
into the appropriately modified Eq. (^9|), in order to 
obtain K^^{uj). 

Approach 3 ("numerical sampling using colored 
noise"): We calculate numerically the time-evolution of 
ps{t) under the action of the stochastic time-dependent 
Hamiltonian defined in Eq. (||), which depends on B(t). 

B{t) is a particular realization of a stationary Gaus- 
sian random process of zero mean and power spectrum 
(BB)^. The density matrix ps{t) has to be averaged over 
a statistical sample of different field configurations B{t). 
This sample is produced by generating the Fourier coeffi- 
cients of B as independent complex Gaussian random 
variables of appropriate variance (given by the power 
spectrum). The field B(t) itself is obtained using a Fast 
Fourier Transform (FFT). After averaging, we may use 
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FIG. 3: The Fourier-transform Kf^{uj) of the equilibrium 
correlator of cyf{t), for different values of the noise strength 
7/(27r) (=0.001, 0.01, 0.05 and 0.1; from topmost to lowest 
graph). The values of the other parameters are: A = 1, As = 
1.2, J = 0.5 and es = 0. Approach 1 and 4: solid line, ap- 
proach 2: dashed line, approach 3: dash-dotted line. 



P{t)P{0) 



(psiidtD), 



(23) 



and Eqs. ^ and (|T|) in order to obtain K^^{lo). To 
this end, the Fourier transform of psii(|t|) is calculated 
numerically, using a FFT on a time-grid of sufficiently 
small step-size Ai and sufficiently large length. The re- 
sults displayed in the figures have been obtained using 
10'' samples and a frequency resolution of Au; = 27r/800. 
The curves have been smoothed by averaging over 5 to 
20 adjacent frequency bins. 

Approach 4 ("weak- coupling equation"): We evaluate 
the analytical results obtained above (see, e.g., Eq. (|2^)) 
with the appropriate numerical values. 

The relevant parameters in our model are es, A5, A, 
the coupling strength J, and the strength 7 of the noise 
force F. We choose the time scale such that A = 1. The 
results discussed in the following have been calculated for 
A5 = 1.2 (5' and B "almost in resonance"). 

To begin our discussion, we note some generic features 
of the results obtained for approaches 1 and 2. Since in 
these cases K^^{uj) is essentially the Fourier transform of 
a density matrix relaxing according to a master equation, 
it consists of several Lorentzian peaks. Their number 
is constrained to be less than the maximum number of 
transition frequencies of the respective system (6 for S + 
B and 1 for S alone, plus possible zero-frequency "pure" 
relaxation). In practice, degeneracies between transition 
frequencies and selection rules reduce that number to 2 
(or 3) for approach 1, and 1 (or 2) for approach 2, for 
es = (or es 0). 

In the limit of weak coupling, J — > 0, all that remains 
is a broadened peak at the transition frequency 2A5 of 
system S alone. In that limit, the results for all three 



models coincide, as expected (see Figs. ^ and||). With 
increasing J, the peaks get broadened and shifted, and 
additional peaks may appear (in the case of approaches 
1,3 and 4). 

Naturally, the behaviour of approach 2 is simplest to 
analyze, since it is the textbook example of a master 
equation applied to a single two-level system. Since the 
correlator (BB)^ is proportional to J^, both the shift 
of the transition frequency and the width of the peak(s) 
increase like J^, for arbitrarily large J. In contrast, the 
dependence of the peak width and the frequency shift on 
the noise strength 7 is non-monotonous. It is determined 
by the evolution of {BB)^ (see Eq. (|o|) and Fig. |) with 
increasing 7. For very small 7, the two-level fluctuator 
B performs very weakly damped oscillations at the fre- 
quency 2 A. Unless it is exactly at resonance with system 
S, the dissipative effects of B on the dynamics of 5' will be 
weak in that regime. The decay rate of S, which is given 
by the power spectrum of B evaluated at 2As, grows 
linearly in 7 (for 7^ < (A| - A2)2/A|). The transition 
frequency of S is shifted upwards or downwards, depend- 
ing on whether the main weight of the spectrum of B is 
located below or above As(A < As or A > As). For 
increasing 7, B performs more strongly damped oscilla- 
tions. The spectrum (BB)^ concentrates around zero 
frequency (see Fig. ||) such that the decay rate of S 
decreases again (like I/7), after having gone through a 
maximum. The magnitude of the energy shift will also 
decrease for increasing 7, simply because the contribu- 
tions of the power spectrum of B lying to either side of 
2As will tend to cancel each other. However, in the limit 
7 — > 00, the shift always saturates at a positive value 
which is independent of A. These facts can be read off 
from the analytical result for approach 2 (written down 
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FIG. 4: The Fourier-transform /('^^(aj) of the equihbrium cor- 
relator of irf (t) , for different values of the coupling strength J 
(=0.1, 0.5, 0.75 and 1.0; from topmost to lowest curve). The 
values of the other parameters are: A = 1, As = 1.2, = 
0.1 and eg = 0. Approach 1 and 4: solid line, approach 2: 
dashed line, approach 3: dash-dotted line. 
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FIG. 5: The Fourier-transform Kf^{u) of the equilibrium cor- 
relator of i3"f (t), for different values of the coupling strength J 
(=0.1, 0.5, 0.75 and 1.0; from topmost to lowest curve). The 
values of the other parameters are: A = 1, As = 1.2, ^ = 
0.1 and es = 1. Approach 1: solid line, approach 2: dashed 
line, approach 3: dash-dotted line, approach 4: dotted line. 



in the special case of es = 0): 

KlAu^) TT-J^ ^' (24) 

s — ±1 ^ ' 

Here the peak width is given by F = 2 j2sffeE(2A5) = 
27r (55)2^^ , the shifted transition frequency is ujq ~ 
2As-2J^'^m'S{2As), and we have defined = {2j+ 
iuj)/{-uj^ + 4A2 + 2iujj). 

The simple master equation (approach 2) is expected 
to come close to the true result as long as the condi- 
tions of the MarkofF and secular approximation are ful- 
filled. This means the coupling strength J has to be 
so small that the resulting decay of S proceeds slowly 
compared with the transition frequency itself (secular ap- 
proximation) and with the correlation time of the bath 
(MarkofF approximation). The latter is given by = 

if 72 < 4A2 and tb = 7^ if 7^ > AA^. 

Now we pass on to a discussion of approach 1, which 
constitutes the exact solution for our physical system. 
The most notable difference to approach 2 is the appear- 
ance of a second peak at the transition frequency 2A 
of the two- level fluctuator B. At small J, the strength 
of this peak grows like J^, while its width is fixed (de- 
pending on 7). In this way, the power spectrum of the 
bath fluctuations shows up in the short-time behaviour 
of the correlator of the system S. This behaviour cannot 
be captured by the master equation (approach 2). In- 
creasing J leads to a frequency shift and a change in the 
width of the "original" peak at 2 As, much like predicted 
by the simpler approach 2. However, in the description 



of approach 1, these changes are due to the change in 
eigenfrequencies and eigenvectors of the combined system 
S + B. At small J, the results of approaches 1 and 2 can 
be shown to coincide using perturbation theory (compare 
also Fig. U and ||, topmost graphs). Deviations from ap- 
proach 2 appear at higher values of J, where the energy 
shift of approach 1 only grows linearly with J (see Figs. 
^ and H, lower graphs). In contrast, the frequency of the 
second peak is suppressed to zero. This behaviour can 
easily be found from the diagonalization of the Hamilto- 
nian for the combined system _B -I- S' in the limit J — > 00, 
when one obtains two pairs of degenerate energy levels, 
separated by 2 J. 

Regarding the dependence on the noise strength 7, the 
same qualitative remarks apply as for approach 2. How- 
ever, it is interesting to note that there is a frequency 
shift with increasing 7 in approach 1 as well (see Fig. 
in spite of the fact that the additional terms in the non- 
secular master equation (|^) seem to describe a purely re- 
laxational dynamics. This is in contrast to the behaviour 
known from the usual form of the master equation, Eq. 
(^) , where the energy shifts can be read off directly from 
the imaginary coefficients in the equation. 

We now turn to approach 3, where the nonlinear bath 
has been replaced by a linear bath (i.e. a colored Gaus- 
sian random process in our case) . Since this model takes 
the full bath spectrum (BB)^ as input, this spectrum 
may also show up in the result for the system correla- 
tor K^^{uj), as is indeed the case. Fig. ^ demonstrates 
that this effect is most pronounced for small values of 7, 
where the bath spectrum has a relatively sharp structure 
(the noise field B{-) acting on S deviates strongly from 
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FIG. 6: The Fourier-transform K'^^{uj) of the equihbrium cor- 
relator of o"f (t) , for different values of the coupling strength J 
(=0.1, 0.3 and 0.5; from topmost to lowest curve). The values 
of the other parameters are: A = 1, As = 1.2, ^ — 0.001 
and es = 0. Approach 1 and 4: solid line, approach 2: dashed 
line, approach 3: dash-dotted line. 

white noise). In these cases, the qualitative agreement 
between approach 3 and approach 1 ("exact solution") 
is much better than that between approach 2 ("simple 
master equation") and approach 1 (see also Fig. Nev- 
ertheless, there are deviations: In particular, there is no 
visible shift of the peaks in approach 3 with increasing 
J. They just become wider and asymmetric (this applies 
especially to the peak at frequency 2 As). For higher val- 
ues of 7, the linear bath (approach 3) in general shows 
less structure than the exact solution, obtained for the 
actual nonlinear bath. 

Finally, we discuss approach 4, where the second peak 
shows up, in contrast to the Markoff approximation. In 
general, we would expect the weak-coupling solution to 
be a bit worse than the simulation of the linear bath with 
colored noise correlations (approach 3), since it is an ap- 
proximation to the latter case. However, the solution for 



the special case eg = turns out to coincide completely 
with the exact solution (approach 1). The solution for 
£5 = 1 (or, more generally, es ^ 0) is good for small 
system-bath coupling J. It fails for increasing J, where 
approach 3 seems to be the better approximation, pro- 
vided es is not too small (see discussion above). 



VI. CONCLUSIONS 

We have discussed a simple model of a nonlinear bath, 
consisting of a single two-level system subject to a classi- 
cal white-noise force. Its action on another two-level sys- 
tem has been analyzed using four different approaches. 
The exact evolution of the density matrix has been com- 
pared with the results obtained using a standard master 
equation, a linear bath (colored noise) which has been 
substituted for the nonlinear bath, and a weak-coupling 
equation. Numerical results for various special cases have 
been discussed. For strong system-bath coupling, the lin- 
ear bath may still be a good approximation in a regime 
where the simple master equation already fails. This ap- 
plies in particular when the bath spectrum has a strongly 
peaked structure. However, deviations between the lin- 
ear and the original nonlinear bath are clearly visible. 

It is hoped that, in a future work, also the (techni- 
cally more involved) case of arbitrary finite temperature 
may be discussed. Further possible extensions include an 
analysis of the higher-order terms in the weak-coupling 
equation as well as replacing i? by a spin of larger mag- 
nitude, in order to observe the transition to the linear 
bath. 
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